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Abstract. The fast solar wind is a collisionless plasma per- 
meated by plasma waves on many different scales. A plasma 
wave represents the natural interplay between the periodic 
changes of the electromagnetic field and the associated co- 
herent motions of the plasma particles. In this paper, a model 
velocity distribution function is derived for a plasma in a sin- 
gle, coherent, large-amplitude wave. This model allows one 
to study the kinetic effects of wave motions on particle dis- 
tributions. They are by in-situ spacecraft measured by count- 
ing, over a certain sampling time, the particles coming from 
various directions and having different energies. We compare 
our results with the measurements by the Helios spacecraft, 
and thus find that by assuming high wave activity we are able 
to explain key observed features of the measured distribu- 
tions within the framework of our model. We also address 
the recent discussions on nonresonant wave-particle interac- 
tions and apparent heating. The applied time-averaging pro- 
cedure leads to an apparent ion temperature anisotropy which 
is connected but not identical to the intrinsic temperature of 
the underlying distribution function. 
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1 Introduction 

It has been known for a long time that the solar wind is a tur- 
bulent plasma with wave activity occurring on a wide range 
of different scales ( Belcher and Davis] 1 1 97 1 1 [Tua nd Marsch, 



1995 Horbury et al. 



2005} . Of course, any plasma wave 
has to fulfill the Vlasov-Maxwell equations, and it can there- 
fore be understood as the space- and time-dependent self- 
consistent interplay between the periodic variations of the 



electromagnetic field and related motions of the particles, be- 
ing represented by their velocity distribution function (VDF). 
In the recent literature, the shaping of distribution func- 



tions due to wave activity has been widely discussed (Wang 
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[etaLl |2006l [Wu and Yoon| [20071 1 Wu et aL| [20091 |Wangand 
|Wu[|2009| . Obviously, the presence of wave forces (or their 
spectra) will lead to a deformation of the distribution func- 
tion with respect to a Maxwellian shape and cause a velocity 
spreading after appropriate averaging over the wave effects. 
The plasma-physics definition of temperature is based on a 
statistical particle ensemble and usually defined in terms of 
the random kinetic energy (via the mean square of the par- 
ticle velocity) in the particles' mean-velocity frame. This 
definition of temperature is not necessarily equal to the ther- 
modynamic temperature, which may be called the intrinsic 
temperature of the particle ensemble. The mean square ve- 
locity fluctuation owing to wave activity is able to cause an 
effective broadening of the distribution function similar to 
real heating, and thus may mimic genuine heating. Therefore 
some authors have referred to this process as apparent heat- 
ing ( |Wang et al.||2006) , others as nonresonant wave— particle 
interactions. The common ground of these wave effects is 
that they are reversible, and therefore not dissipative, and do 
not represent real heating. Collisions, however, might be able 
to dissipate coherent wave motion efficiently, and thus will 
lead to a real heating and an increase in the intrinsic tempera- 
ture ( [Schekochihin et al.| [20081 |Howes||2008) . In the follow- 
ing, collisions are excluded from the treatment of the VDFs, 
in order to demonstrate the collisionless effects of waves and 
reveal the apparent heating due to wave activity. 

Substantial ion temperature anisotropics have been ob- 
served in the solar wind and discussed by different authors 
([Marsch et all [T98T| [20041 |Bale et al.| [20091 |Bourouaine| 
|et al. 20101. Typically the proton temperature is higher per- 
pendicular to the magnetic field than parallel to it. These 
anisotropies have mostly been discussed as being the result 
of the cyclotron-resonant interaction with circularly polar- 
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ized waves, a process which can be quantified by means of the Lorentz force: 



quasi-linear theory (e.g. Akhiezer et al. 1975a Heuer and 



Marsch, 2007 ). On the other hand, such anisotropic distribu- 



tion functions can become unstable if the anisotropy exceeds 



a certain threshold that depends on the plasma beta (Gary 



et al. 2000 2001 ). For a typical solar wind beta of about 1, 



the beta-dependence is not severe and the distribution func- 
tion becomes unstable for T±/T\\ > 2. This instability can 
in turn excite and radiate ion-cyclotron waves. In this way, 
wave excitation can reduce the ion temperature anisotropy 



efficiently and yield moderate and stable values ( Gary 1 1993 
|Bale etal.||2009l . 

The general theory of this wave-particle interaction is well 
established, and in fact many traits of it were confirmed by 
observations (Marsch, 2006). However, the apparent heat- 
ing effects are not well understood and have not yet been 
discussed in the context of ion-wave interactions below and 
near the ion inertial scale. Yet, they should be included in 
an appropriate description of space plasmas that are sub- 
jected to strong Alfven/ion-cyclotron wave activity, such as 
it is typical for the solar wind in the inner heliosphere. To 
study possible physical causes of apparent wave heating is 
the aim of the present paper. First, we discuss the effect of a 
strong plasma wave on an intrinsically Maxwellian distribu- 
tion function. Second, a concise form of the resulting model 
VDF is constructed, and thus we can discuss how such a dis- 
tribution function would look like in a real plasma measure- 
ment made on a space probe. 

2 Model distribution functions 

2.1 Wave effects on the distribution function 

To study the effects that waves have on the shape of a velocity 
distribution function, we have to determine the constants of 
individual particle motion in a given wave field, and then ex- 
ploit the fact that any function of these constants of motion 
is a solution of the Vlasov equation ( [Davidson] |1983| |Stix| 
1992). Here, we discuss the influence of a single monochro- 
matic wave only. It is supposed here to be transversal and 
left-hand circularly polarized, and its magnetic field can be 
assumed to have the form 



(bcos(kz — uit)^ 
bsm(kz-Ljt) 



(1) 



with the constant field component Bq along the z-axis and a 
wave amplitude b. The wave frequency is u and the parallel 
wavenumber k, and its magnetic field is associated with an 
electric field which according to Faraday's law is given by 

1 dB 

elk' 



VxE- 



(2) 



The equation of motion for a single particle with charge q 
and mass m in this electromagnetic field is determined by 



dv 
~dt 



Q_ 
rn 



E H — v x B 

c 



(3) 



It is useful and transparent to write this equation in the com- 
ponents of a cylindrical coordinate system (symmetry around 
the z-axis) as follows: 



dv ± (ui 
dt \k 1 

dt 
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(4) 
(5) 

(6) 



where we made use of the abbreviation <fi — kz — ojt for 
the wave phase. The gyrofrequency is denoted by ft — 
qBo/(mc). A similar set of equations has already been used 
in a test-particle description to describe the nonlinear be- 
havior of particles that are trapped in the wave fields (Mat- 
sumoto 1979|). Also nonresonant heating effects have been 



treated with similar equations under the assumption of low 
plasma betas both in the monochromatic parallel case and 



2006 



for a spectrum of oblique MHD waves ( Hamza et al 
|Lu and Li| [20071 |Li et al.| [2007j |Lu and Chen| |2009| ). How- 
ever, the initial conditions in these cases break the condition 
of the coherent particle motion that is necessary to maintain 
the wave in a self-consistent way. Neglecting the coherence 
and violating the self-consistency may be an appropriate de- 
scription for a minor particle species in the sense of a test- 
particle approach. But the description is insufficient for the 
dominating main species that carry the currents and charges 



maintaining the wave itself. |Li et al. (2007) also consider 



consistent initial conditions that do reflect the coherent wave 
motion and find that these particles do not experience the 
non-resonant heating because they are not picked-up by the 
wave fields. A model for the coherent particle motion of 
the dominating species in a wave field at high plasma be- 
tas should be based on a Vlasov description to take the fi- 
nite thermal width of the distribution into account. Numeri- 
cal self-consistent simulations are another approach to this 



problem (e.g. Li an d Habbal |2005[ |Araneda et al.[ |2008 
|2009| [Maneva et al.[ 2010 ). In order to determine an ade- 
quate model distribution function, we first can determine two 
constants of motion for the kinematic system from Eqs. (|4j- 
[6]). These are the generalized momentum of the particle and 
its total kinetic energy in the wave frame: 



M = v± cos ( 



P = v 2 ± 



UJ 



k (,'n 

2n 



(7) 



(8) 



Both M and P can be shown to be constant, simply by tak- 
ing the derivatives with respect to t and using the equations of 
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motion in 



Eqs. (|4j|6]>. It is known ( Akhiez eretaLl 1 1 975b| ) that 
any distribution function which is a function of these con- 
stants of motion always fulfills the Vlasov equation. There- 
fore, we can make the ansatz 



/ = Nexp 



P 

v 2 " 



exp 



^ 2V V , M 



(9) 



with the normalization factor N and constant coefficients Vth 
and V w , which we can essentially define as the density, ther- 
mal velocity of the VDF and mean fluid-velocity amplitude 
of the particles in association with the wave motion. 

The right second term in the expression for M compen- 
sates the frame-shift of v» by w/k in the definition of P, if 
the condition 

v - = -WT^Jn (10) 

B 1 — w/SZ 

is fulfilled, and if we can ignore the weak effects due to a 
small spread in the parallel direction (fcuii -C f2). Interest- 
ingly enough, this relation then corresponds to the wave po- 
larization relation found by |Sonnerup and Su ( \9G1) for a cir- 
cularly polarized wave with a vanishing parallel bulk drift. In 
their classical solution, they showed that the transversal ve- 
locity is determined by 



oj/k B t 



(11) 



for vanishing drifts in the z -direction. We define the transver- 
sal magnetic field vector as B t — (B x ,B y ,0), with the first 
two components as obtained from Eq. ([T|. The dispersion 
relation of the waves was also given in that work and can be 
written as 

, 9 1 CJ 



whereby the small displacement current in Maxwell's equa- 
tions was neglected. The index j numbers all participating 
species (in the case considered later only protons and elec- 
trons), and lj — c/ujj is the corresponding ion inertial length, 

with the species' plasma frequency uij — A-Krijq 2 /nij. 

We find that after normalization the non-gyrotropic model 
VDF of the particles in response to the wave forces reads 



f w (v±,Vn,cp): 



n 



7r3/2w th 



exp 



V 2 +co 2 /k' 



x exp 



/ 2v ± V w 
! exp R — cos((p — ip) 



(13) 



The first exponential stems from the normalization and does 
not change the structure of the distribution function, but de- 
pends on the particles sloshing velocity amplitude, V^, and 
wave phase speed, ui/k. It is interesting to note, that this 
model distribution function is equal to a Maxwellian distri- 
bution in cylindrical coordinates, yet which is shifted by Vt 



in the transversal direction and can be written in this shifted 
Maxwellian form as 



/~exp 



(V X -V X ) 2 + (Vy-Vy) 2 



(14) 



with the cartesian speed components V x and V y reflecting 
the rigid displacement of the whole VDF in the wave field. 
This is an ansatz commonly used to initialize consistently 
numerical simulations ( |Araneda et al.| |2008| l, and repre- 
sents the sloshing motion of particles in the wave field (e.g. 
|Markovskii etaLl[2009] l. 

A VDF that has an additional intrinsic temperature 
anisotropy is supposed to be represented by a modified bi- 
Maxwellian distribution. The appropriate choice for a VDF 
in a wave field including an intrinsic temperature anisotropy 
is given by 



f a ~ exp 



V 2 



"th-L 



w 2 /k 2 
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thii 
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f2v_ L V Vi 
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(15) 



where different thermal speeds in the perpendicular and 
parallel direction are chosen to account for the intrinsic 
anisotropy. It is important to note, however, that this distribu- 
tion function is not an exact solution of the Vlasov equation 
anymore. Without wave activity, this distribution function 
obtains the usual bi-Maxwellian form. We do not apply this 
modified bi-Maxwellian VDF in the following since we fo- 
cus our considerations on the role of apparent temperature 
anisotropies only. 



(22) 2.2 Wave effects on particle measurements 



The solar wind is permeated by magnetic field fluctuations 
and waves ( |Tu and Marsch| |1995| l of all kind. A particle 
detector that is able to determine the velocity distribution 
function of particles (e.g., such as flown onboard the Helios 
spacecraft) counts particles in different energy and direction 
channels, and thereby integrates the net particle fluxes into 
the various single channels over the so-called sampling time 
T. In a first approximation, we may interpret this time as 
kind of an exposure time (like in photography), and thereby 
neglect time-dependent sampling effects on the instruments 
pointing directions (to different solid looking angles), or on 
the accessibility of the particles to the different energy chan- 
nels during the measurement cycle. Effects of the proper mo- 
tion of the detector with respect to the solar wind flow are 
also neglected for the first estimation. A possibility to han- 
dle the effects arising from this relative motion would be to 
treat the waves as frozen in the solar wind and being with the 
fixed spatial structure convected over the space probe. This 
assumption is called Taylor hypothesis and is valid only for 
& ' ^sw ^ w ' where V$w denotes the solar wind flow ve- 
locity. Since the Alfven speed is typically about a factor of 
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10 less than the flow speed of the solar wind in the spacecraft 
reference frame, this Doppler effect would even increase the 
sampling problem because waves with lower wavenumbers 
appear at higher frequencies for the detector. These lower 
wavenumber structures typically have a higher power than 
waves at higher wavenumbers. Therefore, the relative mo- 
tion of the solar wind with respect to the spacecraft amplifies 
the wave effect on the measured distribution function addi- 
tionally. 

Such a model instrument would not be able to take snap- 
shots of the VDF but integrate the sloshing distribution over 
time T, and would therefore obtain a spread in the VDF due 
to the wave activity. In the following, we determine theoreti- 
cally the influence of this final integration time on the actual 
measurement. 

The relevant spectral range is limited to frequencies that 
are higher than the sampling frequency, 2ir/T, because 
slower motions would be more or less resolved. They would 
merely lead to a rigid shift of the full distribution function 
without deformation. The analysis software of a plasma in- 
strument would set the origin of the reference frame to the 
shifted center of the VDF, so that no change would be de- 
tectable. On the upper frequency side, the acting part of the 
wave spectrum should be limited by the gyrofrequency, be- 
cause the waves are supposed to be strongly damped in the 
case of Alfven-cyclotron waves at this scale, and thus beyond 
it the observed spectral energy goes down significantly. The 
slope of the spectral energy density follows Kolmogorov's 
law at lower frequencies ( |Tu and Marsch|[l995| l. Numerical 
simulations show that at wave numbers around fc«0.8/£ p 
the Alfven/ion-cyclotron spectral slope usually breaks be- 



and at each time. Without any restriction, the spatial posi- 
tion can be taken 2 = 0. Then the time-averaged distribution 
function can be written as 



/ = iVcxp^- 



T 



2v±V w \ sm(mLjT + mip) — sm(mip) 



1 th 



muj 



(18) 



We find that the sum is always zero for all possible positions 
in tp if lo = 2-k/T. So we can simplify the above formula and 
write 



fi 



1 Nexp ( - ^ 11 ) J f 2v± ^ 
\ "th / V "th 



(19) 



and thus we see that the intrinsically non-gyrotropic distri- 
bution function appears to become gyrotropic again after this 
kind of averaging process. A wider spectrum of waves could 
lead to a deformation of the distribution function still in the 
coordinate ip, but as stated above the higher frequencies on 
shorter time scales than the sampling time T are not expected 
to change the result much. 

We assume next an Alfven-cyclotron wave. This wave has 
to fulfill a dispersion relation, which can be taken from the 
cold plasma limit, yielding the dispersion relation as 

(^) 2 -(Hp) 2 + i(H p ) 4 -^(M p )y(fc^)2+4 (20) 



2005 | l. Here, we only assume a monochromatic wave with 
tu = 2ir/T. 

Given all these assumptions, the time -averaging process of 
the VDF under the influence of waves, / w , may be expressed 
mathematically as 



cause of the onset of dissipation at these scales (Ofman et aT] ( |Stix| |1992| |Chandran et al.| |2010[ l. For high values of 

k, the dispersion shows the asymptotical behavior lu — > fl p , 
which can be seen by expanding the function ^/1 + x with 
x = A/(k£ p ) 2 to second order. Typical solar wind parameters 



/wdt. 



(16) 



We can now insert our model function of Eq. ( 13 1, in which 
the last exponential factor can be expanded ( Abramowitz and 



Stegun |1972| l and represented by a sum over modified Bessel 
functions according to the relation: 



acosb . 



: J (a) + 2 I m (a) cos (mb). (17) 



are used for a distance between spacecraft and sun of about 
0.5 AU, which is where strongly non-Maxwellian distribu- 
tions and high wave activity were usually observed in fast 
solar wind by the Helios spacecraft. The plasma beta is set 
to 0.1 and the sampling time to 10 s, which is the real sam- 
pling time of the Helios spacecraft. The constant background 
magnetic field Bq is assumed to be 5 x 10 _4 G, and the pro- 
ton density to be n = 10cm~ 3 . The relative wave amplitude 
is set to b/Bo = 0.25. The distribution function at ip = tt is 
plotted at negative values for v±, in order to make plots that 
can more easily be compared with the observations. The cal- 
culated distribution function is shown in Fig.[T] 



Heavy ions, which are also present in the solar wind (von 



The time dependence of the averaged distribution function / 
is hidden in the wave phase, 4>(t). The coordinates of the 
distribution function (i.e., z, v±, and v») have no time de- 
pendence in this context. The expression for f w in Eq. ( 13 1 
delivers a value of the VDF at each position in phase space 



Steiger 2008 1, react on the wave field in a slightly different 
way. The dependence of V w on the charge-to-mass ratio in 
Eq. ( 10 1 is small for lower frequencies. The thermal speed, 



however, is significantly lower than the proton thermal speed 
by a factor of yjm p /m\. Therefore, the distribution func- 
tion of heavy ions is more narrow than the proton distribu- 
tion function but shifted by the same amount due to the wave 
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Fig. 2. Apparent temperature anisotropy in dependence on the nor- 
malized amplitude b of a wave with uj = 2n/T. The evaluation of 
the anisotropy according to Eq. \2l\ is shown for different values of 
the plasma beta. 



Fig. 1. Distribution function in the presence of a large-amplitude 
wave. The solid line represents the magnetic field direction. The 
velocities are given in units of the thermal speed. The broadening 
in the perpendicular direction is clearly visible. 



motion as can be seen in Eq. (13) . This can lead to a more 
severe deformation of the distribution function after the av- 
eraging process, which might also lead to ring-like apparent 
VDFs for heavy ions in strong wave fields. 

The apparent temperature anisotropy can now be calcu- 
lated by taking the second moment of the distribution func- 
tion according to the formula: 



.4 = 



+ OO+OO _ 

„ / / /iwldw_Ldw|| 

J-± _ -00 
j 1 ,, +00+00 _ 

2 / / fxvf l v±dv±dv\\ 
-00 



(21) 



The integration over ip leads to the factor 2 in the denom- 
inator. The above distribution function of Figure [T] has an 
apparent temperature anisotropy of A — 1.64. 

By varying the wave parameters and the particle thermal 
speed, corresondingly varying forms of the distribution func- 
tion can be created. The dependence of the resulting apparent 
temperature anisotropy on the plasma beta and the wave am- 
plitude b is shown in Fig. [2] 

The obtained distribution function looks from the first 
point of view very similar to a classical bi-Maxwellian dis- 
tribution function of the form 



bm ~ exp 



^th+ 



(22) 








v 



perp 



Fig. 3. Bi-Maxwellian velocity distribution function for a tempera- 
ture anisotropy of 2. 



In Fig. [3] we show a bi-Maxwellian VDF for an (intrin- 
sic) anisotropy of Wtj 1 xA'th|| = 2 ' ^ ne E enera l form of the 
VDF in Fig. [T] can also be approximated by this mathemat- 
ical representation, which underlines the difficulties arising 
from the correct definition of the observed temperature. In- 
terpreting the wave-broadened distribution function / as a 
bi-Maxwellian leads to an apparent difference in the two 
thermal velocities of the bi-Maxwellian VDF. Therefore, the 
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wave-broadening can also be expressed in terms of a corre- 
sponding apparent thermal velocity anisotropy. Applying the 



second moment relation from Eq. (21 1 to a bi-Maxwellian 
distribution function yields the ratio A = v^^/v^u. This 
means that the anisotropy in the apparent thermal speeds is 
given directly by \[A. 

Considering waves with higher frequencies, we must ad- 
mit that the approximation kvii is not valid anymore. 
Then we have to accept further modifications of the VDF in 
the wave field, and we find a non-Maxwellian dependence 
also on the parallel speed coordinate v\\. Accordingly, the 
distribution is 



/h = /wCXp[ — 
u th 



Bo 



' Vw 1 - 



n 



(23) 



If we apply the time-averaging process to this distribution 
function, the results change significantly. The additional 
term makes the distribution function more prolate, already 
at lower wave amplitude and for higher beta values. Thus, it 
becomes even more similar to a bi-Maxwellian distribution 
function as it is shown in Fig. [3] An example for this situ- 
ation is shown in Fig. [4] where we assumed the parameters 
cu = IOtt/T, 6 = 0.1, and f3 = 0.5. 

For comparison, also a typical measurement by the Helios 
2 spacecraft from 1976 is shown. The anisotropy of the ob- 
served distribution function is comparable to the calculated 
apparent anisotropy, whereas other effects such as the formed 
beam along the background field is not reproduced by the 
above calculations. We find that the model distribution func- 
tion fits the observed distribution function well. This shows 
that the broadening effect and the detailed shaping mainly 
depends on the frequency of the waves and the plasma beta. 
The measured distribution function is better represented by 
the corrected distribution function for higher frequencies. 
This distribution is not simply a shifted Maxwellian but has 
a further non-Maxwellian dependence on v\\ , which can rep- 
resent the observations better. 



3 Conclusions 

In this study, we have shown how the VDF is shaped by the 
presence of a large-amplitude wave. In the case of transver- 
sal wave activity, the distribution function obtains a shift in 
the direction perpendicular to the background magnetic field. 
If the distribution function is averaged over time, this shift 
will lead to a smearing in the perpendicular velocity com- 
ponent, which in turn would be interpreted as a temperature 
anisotropy in favour of the perpendicular direction. Every 
real measured distribution function can only be determined 
by sampling within a certain time period, and this implies 
averaging. Thus the resulting temperature as the second mo- 
ment of the VDF reflects this procedure. 

We could demonstrate, using a simplified model, how this 
effect can lead to a significant change in the observed distri- 



bution functions, as plasma measurements are always done 
by counting particles over a certain sampling time T. This 
sampling period corresponds again to time averaging. The 
broadening of particle distributions due to microturbulence 
is a well-known fact, which is exploited in spectroscopy to 
determine remotely, for example, the turbulence level in the 
solar corona (e.g. Kohl et al. 2006) from ultraviolet emission 
line broadenings. In the context of measurements of plasma 
VDF in the solar wind, however, this was not taken into ac- 
count before to the authors' knowledge. 

Also compressive fluctuations can be treated in a similar 
way. However, then broadening would be observed mostly in 
the parallel direction. Consistently with the present empha- 
sis on perpendicular broadening, most recent observations 
show a higher transversal wave activity in almost all cases in 



the fast solar wind (Horbury et al. 2005 1 Alexandrova et al. 
[2008l >. 

The meaning of intrinsic temperature anisotropies should 
be further discussed in the future. As mentioned before, a 
severe limitation to T±/T\\ is observed in the solar wind 
dMarsch et al.l|2Wl[BaieeUll|2009l [Marsch et al.||2009l in 
relation with plasma micro-instabilities which reveal a sensi- 
tive beta-dependence. This finding gives a clear indication 
that the observed temperature anisotropies are largely intrin- 
sic. This finding even more underlines the importance of 
an adequate definition and treatment of the measured plasma 
temperatures. The apparent higher temperature that we find 
in our model VDFs is not the result of resonant heating pro- 
cesses, such as Landau damping or cyclotron resonant wave- 
particle interactions. Hence, this wave-related mechanism is 
reversible, and for a vanishing wave field also the apparent 
anisotropy would disappear. 

Nonresonant wave-particle interactions have been studied 



in the framework of quasi-linear theory by Bourouaine et al. 
(2008). They demonstrate that the nonresonant heating is 
more effective for lower plasma betas since the efficiency of 
the quasi-linear diffusion is mainly proportional to va — v\\. 
We can confirm the beta dependence and find the dependence 
of the parallel particle velocity in the wave frame also in our 
considerations as in Eq. ([8]) for example. However, we could 
show that the dispersion of Alfven/ion-cyclotron waves can 
compensate for the reference frame shift in the Vlasov pic- 



ture. Furthermore, |Bourouaine et al. (2008) found a stronger 
effect of the nonresonant heating on heavy ion species which 
is also consistent with our model. 

In the context of our measurement effect according to 
Eq. dl6}, the determination of temperature is not an ergodic 
measurement anymore if wave fields lead to a coherent par- 
ticle motion as it is the case in the solar wind. This leads 
to the problem that the assumption of Markovian statistics is 
violated since the particle motion is additionally affected by 
the deterministically time-dependent wave motion during the 
averaging. Non-ergodicity in our case means that the tem- 
perature based on time averaging is different from the tem- 
perature based on ensemble averaging. But all real tempera- 
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Fig. 4. left: Like in Fig. [T] but including higher-order corrections, right: Typical proton distribution function measured by the Helios 2 
spacecraft in 1976. The black line indicates the direction of the background magnetic field. 



ture measurements in dilute plasmas have to be done by av- 
eraging over time leading to the apparent deformations in the 
distribution function that we showed above. An appropriate 
ensemble average, however, is not accesible on the required 
scales in the solar wind due to the low particle number den- 
sity. Recently, Hiza nidis et al.| ( |2010[ ) found that the applica- 
bility of classical quasi-linear diffusion is not guaranteed in a 
coherent electromagnetic wave field. They developed a new 
kinetic theory for wave-particle interactions and find time- 
dependent diffusion tensors describing the velocity evolution 
of the VDF. This is a manifestation of the non-ergodicity of 
the measurement process. Maybe a completely different de- 
scription of the microphysical behavior should be applied to 
these coherent cases as it is proposed and discussed in the 
textbook by |Elskens and Escande ( 2003| l. The wave broad- 
ening effect could be excluded locally if, at the position and 
time where the measurement is taken, enough particles are 
present with the local coherent speed additional to their ther- 
mal speed. In all accessible solar wind plasma cases how- 
ever, the number of particles that can be counted under con- 
stant conditions compared to <j> = kz — ojt is practically al- 
ways too small. 

A shorter sampling time would bring the observation 
closer to a "snapshot" of the real distribution function. A 
faster measurement would help avoiding heavy smearing in 
velocity space due to wave activity. The comparatively small 
number density of particles in the solar wind, however, re- 
quires long sampling times or larger geometry factors which 
are not affordable. The Solar Orbiter mission should provide 
new insights, because the improved instruments to be flown 
on board this spacecraft and the higher particle densities ex- 
pected closer to the sun will permit much shorter sampling 
times down to 100 ms. 

In this work, we assumed a monochromatic wave that 
leads to smearing out of the VDF by the time-averaging pro- 



cess. A more realistic assumption would be a broad spectrum 
of waves, which the particles have to follow and the VDF to 
respond to. Furthermore, the particle flux will not be steady 
in all directions during the sampling time T. Also, a realistic 
spacecraft model should be applied, including the detailed 
time dependence of the sampling method as well as the rela- 
tive velocity of the spacecraft with respect to the solar wind. 
If such appropriate models were available, the analysis of the 
distribution function could in turn provide new information 
about the waves, such as their polarization or propagation 
direction. Up to now, in our model only waves propagat- 
ing parallel to the background magnetic field have been dealt 
with. To complete our analysis also oblique wave propaga- 
tion should be taken into account. 
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